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Abstract 



(N 

^ ' The chiral phase transition of QCD is analyzed in a model combining random 

00 . matrix elements of the Dirac operator with specially chosen non-random ones. The 

special form of the latter is motivated by the assumption that the fermionic quasi- 
I zero modes associated with instanton and anti-instanton configurations determine 

the chiral properties of QCD. Our results show that the degree of correlation between 
. these modes plays the decisive role. To reduce the value of the chiral condensate 

\ by more than a factor of 2 about 95 percent of the instantons and anti-instantons 

Q_i' must form so-called molecules. This conclusion agrees with numerical results of the 

^ • Stony Brook group. 



One of the challenging tasks of hadron physics consists in understanding the 
properties of the chiral phase transition in QCD, and its relation to the de- 
confinement phase transition. This fascinating theoretical problem requires 
the consistent combination of several basic concepts of field theory: Dynam- 
ical and spontaneous symmetry breaking, non-trivial topological properties, 
confinement, and a hierarchy of correlation lengths. A deeper understanding 
of the chiral phase transition is also urgently required for the analysis of heavy 
ion collisions at very high energies. Such experiments will be done at RHIC 
and at LHC in the near future. 

Because of its importance, this problem has received considerable theoretical 
attention especially during the last decade. Some papers which are particularly 



relevant for our investigation are listed in Refs. [n]-[14|. We use the notation of 
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Shuryak, Verbaarschot, and Jackson []5|,|7|,[T^ whose papers dominate the field. 
The derivation of the fundamental relations can be found there. 



The starting point for most discussions is the Banks-Casher formula [|T^ 



1/4 



(1) 



It links the value of the chiral condensate, (qq), to the density of the eigenvalues 
/i of the Dirac equation, 

(^i^ + g^^^^^{x)=^^iJ{x), (2) 

at /i = 0. Here, V4 is the four-volume, and A" the gluon field operator. The 
virtuality of typical fermion states is of order Aqcd so that they do not con- 
tribute to p(0). This led to the assumption that p(0) is dominated by very 
special states, the zero modes. Via the Atiyah-Singer index theorem. 
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I d'x f;^f;„ = n+-n^, (3) 



these modes are linked to topologically non-trivial gauge field configurations 
(instantons and anti-instantons). Here, (A^_) is the number of zero modes 
with positive (negative) chirality. The integral on the left hand side contains 
the field tensor and yields the topological charge Q. This charge is -|-1 
(— 1) for an isolated instanton (anti-instanton). To understand the behavior 
of the chiral condensate (and, thus, of hadron masses) as a function of, e.g., 
temperature T it would then suffice to obtain the statistical properties of the 
zero modes associated with QCD-instantons. (For finite instanton separation, 
these are more precisely refered to as quasi-zero modes.) Some time ago, it 
was realized [^,|^ that at T = 0, this statistical problem can be analyzed with 
the help of random matrix models; recently, this analysis was extended for the 
first time to finite temperature ||14||. 



In this paper, we generalize these earlier contributions by investigating the 
dependence of the condensate on two parameters which characterize the tem- 
perature-dependence of the instanton configuration. One of these, is the 
strength of the (diagonal) interaction between fermionic states of opposite chi- 
rality. It is due to the formation of instanton-anti-instanton pairs ( "molecules" ) 
with an optimally attractive relative orientation in color space. This parame- 



ter was also used in Ref. |14[. The second parameter, a, describes the fraction 
of instantons and anti-instantons which form such molecules. The inclusion 
of this new parameter is important because numerical simulations and phe- 
nomenological considerations suggest that it plays a crucial role in chiral sym- 
metry breaking P]. In some of our present work, we actually go beyond this 
simple two-parameter model and investigate distributions of the diagonal ma- 
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trix elements. In this way, we derive some general statements about properties 
of the chiral phase transition, and about the associated critical exponents. 



We see from Eq. that each instanton or anti-instanton is associated with 
one quasi-zero mode. In any given sector of the theory, characterized by the 
value of the topological charge Q = — A^_, the difference between the 
number of instantons and of anti-instantons is given by Q. In this sector, 
the QCD partition function for a theory with Nf fermionic flavors is thus 
approximated by 




qcd(Q)= lim Mjdet i^ + ^?-/° + imJ . (4) 



Here, only quasi- zero modes are taken into account; the matrix in (^) accord- 
ingly has dimension (2A^ + Q). (We define = and assume without loss 
of generality that Q > 0.) The average in (H) extends over all gauge field con- 
figurations A°- with topological charge Q and the usual measure. The symbol 
TUf denotes the mass term. 

In a basis of chirality eigenstates, the operator i^ + g^^"" couples fields of 
different chirality while the mass term in the determinant conserves chirality. 
The determinant can be rearranged so that the first [N + Q) states have 
positive chirality and are localized at the (A^ -|- Q) instanton positions, while 
the remaining N states are the quasi-zero states with negative chirality and 
are localized at the anti-instanton positions. This yields 



The chirality- changing matrix elements have the form 



= <t^l ( + E a^^^^i) + a^KcA 07, • (6) 

Here, we have assumed that the are exact zero modes of the Dirac 

operator when the gauge field A°- in Eq. (j^) is replaced by the field A°'{I]^) 
[A"'{Ii)\ due to a single instanton (anti-instanton). It is assumed that these 
instantons and anti-instantons differ and are well separated, with uncorrelated 
color orientation. The symbol Afluct denotes the fluctuating part of the gauge 
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field not linked to specific instanton configurations. The matrix elements in 
the last line of Eq. (|^) can be replaced by the elements of a random matrix 
with suitable symmetry properties PJ7|,p^. This leads to a random- matrix 
model with partition function 



Zp{Q) = J V[W] exp 



-^^tr(iyiyt) 



TV 



X 



/=i 



ndeti, , '^'::'^). (7) 



The matrix TV is a rectangular matrix with (A^ + Q) rows and columns. 
Depending on the underlying symmetry, we distinguish three random matrix 
ensembles pertaining to the partition function labeled by /3 = 1, 2, and 
4: the Gaussian orthogonal ensemble (GOE, (3=1), the Gaussian unitary 
ensemble (GUE, /? = 2), and the Gaussian symplectic ensemble (GSE, /3 = 4) 
with real, complex, or quaternion real matrices W, respectively. As shown by 
Verbaarschot 0, these three cases correspond to the following gauge groups: 
SU(2) in the fundamental representation for /? = 1, SU(A"c) with A^^ > 3 
in the fundamental representation for (3 = 2, and non-Abelian gauge groups 
SU(A'c) for all Nc with fermions in the adjoint representation for (3 = 4. 
In this letter, we are concerned with the case (3 = 2 only. The integration 
measure in Eq. (|^) is the Haar measure. The parameter S is a real number 
which is related to the value of the condensate at T = and which we take 



to be independent of temperature. Lattice calculations ||T^ suggest that S 
can depend on temperature only weakly; such a dependence cannot affect our 
results quahtatively. 

The matrix D is real and diagonal. As mentioned above, D is included to 
describe the non-random components of the Dirac matrix elements due to 
the formation of instanton-anti-instanton pairs with an optimally attractive 
relative orientation in color space. The matrix D is diagonal because only 
the two quasi-zero modes of opposite chirality belonging to such a pair are 
coupled by a non-random matrix element. We expect the elements of D to grow 
with increasing temperature. This point will be substantiated below. The T- 
dependence of D generates the T-dependence of the chiral condensate. In the 
sequel, we study several distributions of the matrix elements of D. Analytical 
results are obtained only for the simplest of these, for which a fraction a of the 
elements of D is non-zero, all non-zero elements being equal. We believe that 
the distributions we study should account for the most important features 
of QCD related to the chiral phase transition, with one possible proviso. It 
would seem much more realistic to replace the random matrix by a banded 
random matrix [Q. This would account for the fact that the interaction of 
widely separated zero modes should be suppressed, and would introduce a 
correlation length in the theory. Work on this problem is under way. 
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We turn to the details of the calculation. The determinant in can be 
rewritten as an integral over ant i- commuting (Grassmann) variables [|l^, 



Z2{Q) = V[W]V[ij]exp 



Nf 



-Ni:hT{ww^) 



The ( ant i- commuting) vector ip has Nj {2N + Q) components which we label 
as follows. The upper index / = 1 . . . Nf denotes the flavor, and the lower index 
numbers the (anti-) instanton zero modes. The first N + Q components of the 
lower index denote the instanton zero modes labeled by Ik with k = l . . . N+Q, 
and the remaining components the anti-instanton zero modes labeled by /; 
with / = 1 . . . A^. Performing the integration over W we obtain 



Z2iQ) 



v\ 



exp 



(9) 



Summation over repeated indices is assumed, and we have set Du = Di. The 
four-fermion terms can be removed by means of a Hubbard-Stratonovitch 
transformation at the expense of introducing 2Nj additional integration vari- 
ables. These can be taken to be the elements of a complex Nj x Nf matrix 5*. 
The result can be rearranged as 



Z2{Q) 



vism 



exp 



N+Q 

k=N+l 
'S + M 



N 



k=l 



,(10) 



where the T/'-vectors are now in flavor space, the Dk are multiplied by the 
Nf X Nf unit matrix, and the matrix M is a diagonal matrix with entries rrif. 
The integration over the Grassmann variables can now be performed to yield 



Z2{Q) = J V[S] exp [-Ars2tr(^5^)] det^{S + M) 



N 

X W det 

fc=i 



S + M 



5 



The (gg)-condensate is given by 
1 d loe: Z 



ViNf DM 



M=0 



which, using ([TTp, can be expressed as 



'12) 



N 



V[S] exp [-ArE2tr(^^t)j det^(^) [] det .^^ '^^^ 



4iV/ 



Qtr(5 



-In 



k=l 



(13) 



In the hmit N ^ oo, the integral in (|TT|) can be evaluated in saddle-point 
approximation. Neglecting terms of order Q, we obtain 



TV 



k=l 



NJ:^S = J2iS + M) {S + M){S^ + M) + D 



-1 



(14) 



Any square complex matrix of dimension Nf can be written as 5* = UAV ^, 
where A is a diagonal matrix with real and non-negative entries, U is unitary, 
and Ve\JiNf)/(^Zi U(l). The solution of (|1§ in the limit of zero masses is, 
thus, obtained by solving 



N 



1 



-y 



(15) 



for A. The (gg)-condensate then follows from Eq. 
1/4 

where njnst is the instanton-anti-instanton density. 



S^A — ^ , ^ — S^Aninst ; 



(16) 



In a realistic physical situation, the entries Dk are distributed according to 
some (normalized) distribution P{Dk). In the limit — > 00, we can map the 
(discrete) interval = 1 . . . A^ onto the (continuous) interval x = . . . 1. The 
desired value of A is then obtained by solving the equation 



dx 



A2 + D^x) 



(17) 



for A. Here, D{x) is defined by 

D{x) 

J P{y) dy = x. 
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Although the above equations cannot in general be solved analytically one can 
still make interesting statements about the connection between the existence 
of a condensate and the analytical form of P{y) for small y. A condensate exists 
if Eq. (p!?!) has a positive solution for A. This is the case if Jq dx/D'^{x) > S^. 
Since D{x) increases monotonically by definition this condition is automat- 
ically fulfilled if D{x) ~ x"' for small x with a > 1/2. Using Eq. (|18D, this 
is equivalent to P{y) ~ y'^ for small y with b < 1. The distribution P is a 
function of temperature. Hence, necessary conditions for the vanishing of the 
condensate above some critical temperature are 



Piv) ~ for y -* with 6 > 1 
1 



da; 
DHx] 



(19) 
(20) 



While these conditions set non-trivial constraints on the general distribution 
of the diagonal elements, it is more convenient for our purposes to consider a 
number of interesting special cases where analytical results can be obtained. 
The simplest possible case is that where all elements of D are equal, corre- 
sponding to P{y) = 6{y — d) and D{x) = d. The solution for A can be obtained 
immediately to yield 



{qq) = SninstVl " i^d^ . (21) 
This version of the model with d = ttT was considered by Jackson and Ver- 
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As a more realistic case we assume that a fraction a G [0, 1] of the elements of 
D are equal to a common value d and the remaining elements are zero. This 
corresponds to P{y) = o.5{y — d) + (1 — a)6{y), D{x) = for < a; < 1 — a, 
and D{x) = d for 1 — « < x < 1. Again, we can solve for A to obtain 

{qq) = ^^1 - (Srf)2 + .^[1 - (Srf)2]2 + 4(Srf)2(l - a) . (22) 

This result is shown graphically in Fig. 1. Its most important property is that 
(qq) can approach zero only if a is very close to one. This observation is in strik- 
ing agreement with the results of a specific model by Ilgenfritz and Shuryak |^ 
and with the general observation that the chiral condensate changes abruptly 



at the chiral phase transition |11] while the instanton density changes only 
very smoothly. Note also that the relatively sharp drop of the chiral conden- 
sate close to the critical temperature is in qualitative agreement with results 
obtained by the Bern Group in the framework of chiral perturbation theory 



I TBI. We conclude from our results that this is a generic property which follows 



7 



from the symmetries of the interaction and the fundamental properties of in- 
stanton molecules. Thus, any specific model should reproduce these features. 




Fig. 1. Plot of the {qq) condensate in units of its value at T = as a function of d 
and a defined in the text. Both d and a are functions of temperature so that the 
dependence of {qq) on T corresponds to a specific path on this surface. 



At this point, we are in a position to discuss specific numbers corresponding to 
the real physical situation. The QCD vacuum state is characterized by a num- 
ber of quantities, one of them being the gluon condensate (fG^^G^^), which 
in turn can be related to the instanton-anti-instanton density (in Euclidean 
space-time). Phenomenologically, ninst ~ 1 fm~^ at T = 0. It is expected 
that this value changes only little with temperature up to Tc. The typical in- 
stanton radius p is equal to the color correlation length which is determined 
by many phenomenological considerations (e.g., from pomeron phenomenol- 
ogy or from the gluon form factor |T^) to be about 0.3 fm at T = 0. The 
same number is determined by fitting a host of hadronic properties in the in- 
stanton vacuum model of Shuryak and co-workers. These numbers are further 
corroborated by specialized lattice studies using cooling techniques to single 



out the quasi-classical field configurations pO| , pT| . Furthermore, the separation 
of instanton and anti-instanton in a correlated pair is i? ~ 1 fm |20], and the 
empirical value of the quark condensate is {qq) ~ —(225 MeV)^ 



The diagonal matrix elements d have been calculated by Diakonov and Petrov 
P|, Eq. (22). For an optimally attractive orientation of instanton and anti- 
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instanton this implies 



\d\ = -^^ 0.4 fm"^ (23) 

at T = 0. On the other hand, we have for the ground state with no noticeable 
instanton-anti-instanton correlation (a — »• 0) 



^ -1.5 fm . (24) 

T=0 



These numbers imply y 1 — (Srf)^ ~ 0.8 at T = 0. Thus, it is not sufficient 
that a approaches one, but in addition d has to increase by about fifty percent 
as T increases from zero to Tc. This can be achieved by an increase of p and/or 
a decrease of R. Lattice studies indicate that p is approximately independent 



of T ||21[]. Recent calculations by T. Schafer and E. Shuryak in the framework 
of the instanton liquid model show that R decreases with T Thus, an 
overall inrease of d with temperature is plausible. 



We can compute various critical exponents from the above results. Our order 
parameter is (qq), and the control parameter is the temperature T. Let us 
concentrate on Eq. (|2^). The two parameters d and a depend on T in a way 
which is beyond the scope of our model. However, some qualitative statements 
can be made on physical grounds. At T = 0, d has some finite value whereas 
a is zero. Both parameters will grow with temperature, but a is bounded by 
one. It is reasonable to assume that d is an analytic function of T, but no such 
assumption can be made for the behavior of a(T) close to a = 1. There is even 
the possibility that a changes discontinuously to one resulting in a first order 
transition. We note that our model is rich enough to allow for that possibility 
but will assume in the following discussion that a is a continuous though not 
necessarily analytic function of T. 

It is clear from Eq. ( |22D and Fig. 1 that necessary conditions for a vanishing 
condensate are Ed > 1 and a = 1. Let us denote the temperatures at which 
these two conditions are first satisfied by and Tq,, respectively. The larger 
of these two temperatures is the critical temperature Tc, but we do not know 
a priori which one it will be. In principle, there are three possibilities: (i) 
Tc = Td > Tq,, (ii) Tc = Td = Tq,, and (iii) Tc = Tq, > T^. We are interested in 
the critical exponent jS defined by 

{qq) ~ \tf with t = ^ . (25) 

J- c 

(This j3 is, of course, different from the (3 labeling the random matrix ensem- 
bles.) For cases (ii) and (iii), jd depends on the specific form of «(T) close to 
Tq. We parameterize this dependence by 1 — a{t) ~ [tY' for t — > 0^. We then 
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obtain for case (i) [3 = |, for case (ii) j3 = \ ioi x > 2 and /9 = | for x < 2, 
and for case (iii) /? = f • 

The critical exponents 7 and 5 are defined by 



1+1-7 -fv, 
X ~ |t| with X ~ 



(26) 

m=0 



dm 

{qq) ~ m^/^ at T = T, . (27) 

In this case, the mass m plays the role of the conjugate field breaking chiral 
symmetry explicitly (the analog of the external magnetic field in the case of a 
ferromagnet). Both 7 and 6 can be computed in a straightforward manner by 
going back to Eq. (p^. The results, however, again depend on the relationship 
between and T^. We find that in case (i) 7 = 1 both above and below the 
transition and 5 = 3 in agreement with Ref. |T^. In case (iii) 7 = both 
above and below the transition and 6 = 1. Case (ii) is more complicated. We 
obtain that above the transition 7 = 1 whereas below the transition 7 = 1 for 
X > 2 and 7 = | for x < 2. Furthermore, 6 = 3 here. In all three cases, the 
Widom scaling law, 7 = /?(5 — 1), is obeyed. 

Case (i) and case (ii) for x > 2 yield the standard mean-field exponents while 
case (ii) for a; < 2 and case (iii) give rise to critical exponents which do 
not seem to be in a familiar universality class. This finding deserves further 
discussion. 

When considering a second order phase transition, one usually assumes that 
the "input" parameters of the system (in our case, d and a) depend on the 
control parameter in a smooth and analytic way. The non-analytic features of 
the "output" parameters (e.g., the order parameter) at the critical tempera- 
ture are a consequence of the phase transition. It is not in the spirit of such 
considerations to introduce non-analyticities at the outset, i.e., in the input 
parameters. In our specific case, we are not able to specify how a approaches 
one as T ^ T~ . However, a is special in the sense that it is strictly equal to 
one for T > T^. Thus, requiring a{T) to be continuous and differentiable at 
Tq, yields 1 — aiT) ~ (Tq — T)^ for T — > T~ and, hence, x = 2 in the above 
results for (3 and 7. Therefore, if it should turn out that either Td > or 
Td = Ta with a(T) continuous and differentiable at T^,, we reproduce the mean 
field exponents of a second order phase transition. This would be consistent 
with the recent discussion in p4[. If Td < or Td = T^ with a{T) continuous 



but not differentiable at we obtain a second order phase transition with 
non-standard critical exponents. Finally, if a is discontinuous at Tq, > T^, we 
obtain a first order transition. 

In summary, we have extended the chiral random-matrix model of Refs. |5|,[^,p!^ 
in such a way that the model allows either for a first-order or for a second-order 
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chiral phase transition. For the second-order transition, we have calculated the 
critical exponents, and we found a number of possible solutions, all consistent 
with universal scaling laws. A decision between them is beyond the present 
scope of our model. We believe, however, that these options define the generic 
possibilities of any physically realistic model of the chiral phase transition. 
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